Real-Space Renormalization Group Study of the Two-dimensional Blume-Capel 

Model with a Random Crystal FieldB 
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The phase-diagram of the two-dimensional Blume-Capel model with a random crystal field is in- 
vestigated within the framework of a real-space renormalization group approximation. Our results 
suggest that, for any amount of randomness, the model exhibits a line of Ising-like continuous tran- 
sitions, as in the pure model, but no first-order transition. At zero temperature the transition is also 
continuous, but not in the same universality class as the Ising model. In this limit, the attractor 
(in the renormalization group sense) is the percolation fixed point of the site diluted spin-1/2 Ising 
model. The results we found are in qualitative agreement with general predictions made by Berker 
and Hui on the critical behaviour of random models. 
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I. INTRODUCTION 

The Blume-Capel (BC) model is a spin-1 Ising model, 
originally proposed to study first-order magnetic phase 
transitions jlj: its phase-diagram presents a line of con- 
tinuous transitions and a line of first-order transitions, 
separated by a tricritical point. The Hamiltonian of the 
model is given by 



Hbc = -J S - S - ' A X S ' 



(1) 



<hj> 



where the first sum is over all nearest-neighbor pairs on 
a lattice and the last one is over all sites, J is the ex- 
change constant, A is the crystal field and Si = ±1,0. 
Later, a generalization of the BC model was introduced, 
the Blume-Emery- Griffiths (BEG) model ||]: it has been 
used to study a rich variety of physical systems, among 
them 3 He- 4 He mixtures. Its Hamiltonian reads: H = 
Hbc — KJ2<i j> SfSj and the parameters J, K and A 
were originally related to the energy interactions between 
the constituents of the system. In 3 He- 4 He mixtures, the 
state S = represents a 3 He atom while 4 He atoms are 
denoted by S = ±1 states: the superfluid transition cor- 
responds to the symetry breaking between the ±1 states. 

More recently, the critical behavior of 3 He- 4 He mix- 
tures in random media (more precisely, in aerogel) has 
been modeled by a BC model with a random crystal field 
(RFBC). The presence of the porous media is taken into 
account by the introduction of a site-dependent crystal 
field, which follows the probability distribution 



P(Ai) = P S(Ai - AO + (1 - P ) 5(Ai - A 2 ), 



(2) 



where Ai is the field at the pore-grain interface and A 2 is 
a bulk field which controls the concentration of 3 He atoms 
. More precisely, the BEG model has been used to de- 
scribe 3 He- 4 He mixtures, and the biquadratic exchange 
parameter, K, is related to the interaction energy be- 
tween a He-^He atoms, V a p, through V33 + V44 — 2V34. 
Since V a p is nearly independent of a and /3, one can as- 
sume that K is zero, regaining the Blume-Capel model. 

From the theoretical point of view, the presence of 
randomness may affect the critical behavior of systems 
in a drastic way. Random bonds || and random fields 
Q effects on phase transitions have been studied for a 
long time. Briefly, the effect of random fields on multi- 
critical phase diagrams is the following: the presence of 
an infinitesimal ammount of randomness eliminates non- 
symmetry-breaking first-order transitions and replaces 
symmetry-breaking first-order transitions by continous 
ones in two dimensions (d = 2), while for d > 2 tricrit- 
ical points and critical end points are depressed in tem- 
perature and first-order phase transitions are supressed 
only at a finite amount of the disorder . Whether the 
first-order transition in two dimensions is replaced by a 
continuous transition in the same universality class as 
the spin-1/2 Ising model, or otherwise, is still an open 
problem [Q. 

To the best of our knowledge, theoretical works with 
the RFBC model have used some sort of mean-field-likc 
approximation ||[|||JI| . These approximations describe 
correctly the behavior of high-dimensional systems and, 
even in those cases, a meaningful discussion on universal- 
ity classes is not possible. The so called effective-field ap- 
proximation p|, for instance, cannot describe first-order 
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phase transitions; on the other hand, standard mean-field 
approaches (which assume that each spin interacts with 
all other spins in the system) do not lead to a correct 
discussion of percolation effects. 

In the present work, we employ a real-space renor- 
malization group (RSRG) approximation to discuss the 
two-dimensional RFBC model. Our approximation takes 
into account spin correlations at all levels, and allows for 



the discussion of first-order transitions |10|, universality 
classes, multicritical points |11|, etc. The crystal field 
probability distribution we chose is slightly different from 
Eq. (0), namely 



V(Ai) = p S(Ai + A) + (1 - p) 8(Ai - A). 



(3) 



We believe that the important physical ingredient lies on 
the presence of randomness and not on the exact form 
of the probability distribution. In fact, we performed a 
mean-field calculation using Eq. (|§|) as the probability 
distribution and the results are qualitatively equivalent 
to those obtained using different distributions and the 
same mean-field approximation. 

The remainder of this paper is organized as follows. 
In section II we outline the formalism and discuss some 
technical points, in section III we present the results, and 
in the last section we summarize our main conclusions. 



II. FORMALISM 

We approximate the Bravais lattice, which in our case, 
is the square lattice, by an appropriate hierarchical lat- 
tice. We chose one of the simplest cells, depicted in Fig. [|; 
albeit its simplicity, this cell has been used with success in 
the study of many ferromagnetic systems. We note that 
the results obtained are exact on the chosen hierarchical 
lattice but only approximate on the square lattice. In 
particular, one does not expect to obtain results as pre- 
cise as those from Monte Carlo simulations or conformal 
invariance arguments. Nevertheless, universality classes 
and the order of the transitions are very well described by 
RSRG approximations, particularly in two dimensions. 




J',K' 



FIG. 1. Construction of a hierarchical lattice adequate to 
simulate the square lattice. Si and S2 denote terminal spins 
while S3 and S4 denote internal spins. The original lattice is 
depicted on the left-hand side of the figure, with parameters 
p, J, K, and A; after doing a partial trace over spins S3 and 
S4 we are left with the "new" lattice (right-hand side), with 
renormalized parameters p' , J' ,K' and A'. Note that the con- 
struction of the hierarchical lattice is done by reversing the 
length scale transformation. 



We then impose that the correlation function between 
the two terminal sites of the original and renormalized 
graphs are preserved 113] : 



exp(-/37Yi 2 ) = Tr exp(-/?Wi i2 ,3,4), 



(4) 



where Tr means a partial trace over the internal sites 
of the cell (S3 and £4 in Fig. |l|). We rewrite the cell 
Hamiltonians as a sum of "bond" terms (from now on, the 
factor P will be absorbed into the interaction parameters) 



Hi. 



2,3,4 



'J(SiSs + S1S4 + S3S4 + S2S3 + S2S4) 
K{SlSl + SjS; 



Ai 



q2 q2 

A 4 



q2 q2 
°2 °3 



S 2 S 4 ) 



A 4 



+ (2^ + 3=^ + 3=^| + 2=^), (5) 

where the crystal field Aj follows the probability distri- 
bution given by Eq. (||) , and 

W12 = -J'StS 2 - K'SfSl + (^-Sf + =^Sf ) + G", 

(6) 

where primed quantities are renormalized parameters 
and G' is a constant, generated by the renormalization 
procedure. We comment below on the renormalized prob- 
ability distributions. 

Note that this way to write the cell Hamiltonians (Eqs. 
H and ||) is equivalent to attribute weights to the sites 
in the one-site (crystal-field A) interaction, according to 
their coordination number. This is necessary for finite 
lattices in order to approximate correctly the infinite lat- 
tice behavior (see, for instance, Ref Jl^]). By using the 
above procedure we obtain the exact value for the point 
where both ferromagnetic and paramagnetic phases co- 
exist at zero temperature for the pure (p = 0) Blumc- 
Capel model on the square lattice, namely (A/J) c = 2 
(see Fig. |). 
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FIG. 2. Phase diagram of the pure Blume-Capel model, 
where k is the Boltzmann constant and T is the temperature. 
The first order (dashed) line flows to a zero temperature fixed 
point, where the largest eigenvalues for odd and even sectors 
of the RGT matrix equals b d , where 6 is the length-scaling 
parameter and d is the dimension of the system. TC is the 
tricritical point, which flows to J* = 1.51, K* = 0.051 and 
A* = 3.01 (within the present approximation). The contin- 
uous line, to the left of TC, is atracted to the spin- 1/2 Ising 
fixed point. O (D) stands for ordered (disordered) phase. 

Some points are worthy stressing at this stage. First, 
we comment on the presence of the biquadratic inter- 
action K in our formalism. Although we are treat- 
ing the Blume-Capel model (K = 0), the parameter 
K is generated by the renormalization procedure and 
it must be taken into account to follow the renormal- 
ization path. To restrict oneself to a subspace which 
is not invariant usually leads to spurious results. Sec- 
ond, the renormalization procedure will introduce ran- 
domness in all renormalized quantities {J',K' and A'). 
One possible approach would be to follow the successive 
renormalized distributions of these parameters in order 
to study the phase diagram. We adopted an alternative 
way which forces the renormalized distributions to be the 
same as the initial ones, but with renormalized parame- 
ters, namely, V' ap {J) = S(J - J'), V' ap (K) = S(K - K 1 ) 
and V' ap (Ai) = p' S(A t + A') + (1 - p') 5(A< - A'). The 
values of J' and K 1 are obtained by imposing that the 
first moment of the actual distributions for J and K and 
of V ap (J) and T" ap (K) are equal, respectively. The val- 
ues p' and A' are calculated imposing that the two lowest 
moments of V' ap (A) match those of the real distribution. 
This procedure has to be used with some care: in some 
systems where the random-field mechanism is important 
and the initial randomness is on the interaction (J, for in- 
stance) , forcing the field back into a uniform distribution 
leads to incorrect results. In Ref. Q|, for instance, the 
crystal-field probability distribution is maintained uni- 
form throughout the renormalization procedure. Conse- 
quently, the random model critical behavior is character- 



istic of a high-dimensional system: the critical temper- 
ature of the tricritical point diminishes as randomness 
is increased but only reaches the zero temperature axis 
at a finite value of the disorder. As thouroughly dis- 
cussed in Ref. , the mechanism responsible for the lack 
of first-order phase transitions in two-dimensional ran- 
dom systems is the disorder in the crystal- field, which 
is not taken into account by approximations such as the 
one used in Ref. 113]. In the model we study in this 
paper, however, the important physical ingredient is the 
disorder in the field, which is not approximated away by 
our RSRG procedure. Finally, we would like to mention 
that the way we treated the random field distribution is 
not unique. In this work we assume that only one field 
acts on each site and a weight is associated to the fields 
(this weight is the ratio between the coordination num- 
ber of the site in the cell and the coordination number 
of the site on the square lattice). Conversely, one could 
also take the number of fields acting in a given site as 
equal to the coordination number of the site in the cell. 
We performed calculations using both procedures above. 
The results do not vary qualitatively (and some times 
quantitatively) from one approach to the other. The ap- 
proach we chose, however, leads to simpler expressions, 
which are easier to deal with analitically. 

The expressions connecting renormalized and original 
parameters are easily obtained following the procedure 
outlined above but are too lenghty to be explicitly written 
here. Formally, they can be expressed as 

f = J'(p, J, K, A); K' = K'(p, J, K, A); 

A' = A'(p, J, K, A); p' = p'(p, J, K, A). (7) 

Critical points are then evaluated as non-trivial fixed 
points of the above relations; phases are identified ac- 
cording to the attractor of their points. The order of the 
transition is obtained through the study of the largest 
eigenvalue of the rcnormalization-group transformation 
(RGT) matrix 0. More precisely, a first-order phase 
transition such that m =< S > is discontinuous at the 
transition point is signaled by the presence of an eigen- 
value equals to b d in the odd sector of the linearized RGT 
matrix, where b is the length-scaling parameter and d is 
the dimension of the system. If q =< S 2 > is discontin- 
uous, the b d eigenvalue belongs to the even sector of the 
RGT matrix. In the present case, b d — N/N' — 5, where 
N is the number of bonds of the original cell and N' is 
the number of bonds of the renormalized one. 

III. RESULTS 

In Fig. | the pure (p = 0) phase diagram is depicted, 
for completeness. We would like to stress that the dashed 
line (and its zero temperature point) is atracted to a fixed 
point (which depends on the approximation) where the 
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largest eigenvalue for both even and odd sections of the 
RGT matrix equals to b d , indicating a first order phase 
transition in m and in q. Note that the K = plane is 
not an invariant one and the biquadratic interaction K 
is generated by the renormalization transformation. 

Following Ref. , the first-order transition should van- 
ish for p > (random model). This is actually the be- 
havior we observe. In fact, the first-order fixed point 
attractor of the dashed line in Fig. ^ is found to be un- 
stable along the p direction. This is the expected physical 
behavior when randomness is introduced. On the other 
hand, the attractor of the pure Ising model transition 
line, namely (p* = 0, J* = 0.4407, if* = -0.0731, A* = 
-co), is stable along the same direction. There are still 
two possibilities for the random model critical behavior: 
either the whole line of continuous transition belongs to 
the universality class of the spin- 1/2 Ising model or an 
unstable fixed point at finite temperature separates the 
Ising critical line from another continuous line which be- 
longs to a new universality class. Our results suport the 
first option: the Ising critical line extends down to the 
zero temperature point (see Fig. [|, where typical phase 
diagrams for p < p c — 1/2 are depicted). Here there are 
still two possible scenarios. The continuous transition for 
p or 1 belongs either to the pure or to the disordered 
Ising model universality class. For the hierarchical lattice 
we use in this work, the specific heat critical exponent of 
the pure Ising model, a, is negative and disorder is ir- 
relevant, according to the Harris criterion [g. Therefore, 
the continuous transitions depicted in Figs. 3, 4 and 5 
belong to the pure Ising model universality class. For 
the corresponding model on a two-dimensional Bravais 
lattice, where a = 0, the Harris criterion is inconclu- 
sive. The accepted behavior, when disorder is present, 
is the following: critical exponents of the random model 
retain the same values as their pure conterparts but log- 
arithmic corrections are introduced by randomness [|15| . 
Experimental results also indicate the same critical ex- 
ponents for the pure and random two-dimensional Ising 
model fTq| . On the other hand, when a is positive, as 
in the three-dimensional Ising model, disorder makes the 
system to crossover to a new universality class. 

Note that the critical value of A/ J which separates 
the ordered and disordered phases, (A/J) c , increases as 
p decreases, in contrast to the result obtained by a cluster 
variational approach on a similar model [Q , which leads 
to a constant value for (A/J) c for any p < p c . The 
latter result might be an artifact of the cluster variational 
approximation. 
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FIG. 3. Phase diagrams of the RFBC model for p < p c (see 
text), (a) stands for p — 0.1 and (b) stands for p = 0.3. There 
is no first order transition and both continuous lines belong to 
the pure Ising model universality class. The frontier at zero 
temperature is attracted to the random fixed point (see text). 
O (D) stands for ordered (disordered) phase. 



At zero temperature, points on the frontier between 
the disordered and ordered phases flow to a random fixed 
point, (p* = 1/2, J* = oo,K* = |ln(2) - J*, A* = oo), 
such that J* /A* = 0. This is the percolation fixed 
point of the site-diluted spin-1/2 Ising model. In fact, 
for A = oo the RFBC model is equivalent to the random 
site spin-1/2 Ising model, where sites are present or ab- 
sent with probability p or 1 — p respectively. This comes 
from the fact that, for A = oo, a +A crystal field acting 
on a given site forces that site to be in the S = state 
(absent), while a — A field forces the site to be either in 
the state S = 1 or in the state S = — 1 (both represent 
a present site). Thus, only for high enough p an infinite 
cluster of S = ±1 states will form and will be able to 
sustain order. Exactly at p — p c , there is such an infinite 
cluster but its critical temperature is zero. Therefore, 
the critical parameter (A/J) c reaches oo for p = p c (see 
Fig. [|). Our evaluation of p c is 1/2, while the accepted 
value for the site percolation critical probability on the 
square lattice is p c — 0.5927 jL7j. It is not unusual that 
small-cell RSRG approximations fail to obtain a quan- 
titatively precise value. Note, however, that we do ob- 
tain the correct qualitative behavior, i.e., a finite value of 
p c (contrarily to the standard mean-field approximation, 
which predicts p c = ||). 
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FIG. 4. Phase diagrams of the RFBC model for p = 1/2, 
which is the value of p c in our approximation. The critical 
line touches the zero temperature axis at A/ J = oo. O (D) 
stands for ordered (disordered) phase. 

For p > p c , the critical line never touches the (A/ J) 
axis. Even at A/ J = oo the infinite cluster of S = ±1 
spins is, on a large scale, a two-dimensional object and 
its critical temperature is finite (see Fig. ||). 




A 
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10 20 30 

FIG. 5. Phase diagram of the RFBC model for 
p = 0.6 > p c - The critical line never touches the zero temper- 
ature axis. O (D) stands for ordered (disordered) phase. 

At this point, it is worthwhile to compare our re- 
sults with those from mean-field calculations (see Refs. 
HHH). Standard mean-field analysis leads to a first- 
order transition inside the ordered phase, ending in a crit- 
ical end point, a reentrant behavior in the kT/J x A/ J 
diagram and a physically incorrect value for p c . We have 
already commented on this last feature. Concerning the 
first-order transition inside the ordered phase, it has been 
shown that it is unstable against randomness in two di- 
mensions 0. Thus, it is expected that a reliable ap- 
proximation to a two-dimensional system will not find 
such transition. Finally, we found no reentrance in our 
results; actually, in some other models reentrant behav- 
ior has been found for d — 3 systems, but not in their 
two dimensional counterpart (see, for example, Ref. psf ). 
We should also point out that more sophisticated mean- 
field-like approximations have been applied to the RFBC 
model. They lead to a finite value of p c but still pre- 
dict the existence of first-order transition in the random 
model as well as a reentrant behavior. Hence, the results 
shown in this work reflect the correct qualitative behavior 
of the RFBC model in two dimensions. 

IV. SUMMARY 

A RSRG procedure is applied to the RFBC model in 
two dimensions. Our calculation recovers the correct 
phase diagram of the pure model and predicts that no 



first order phase transition is maintained when random- 
ness is introduced. This is in accordance with general 
predictions for two-dimensional disordered models [Q. 
We also obtain that the whole line of continuous tran- 
sitions, for j) ^ 0, belongs to the Ising universality class, 
discarding the existence of an unstable fixed point at fi- 
nite temperature. The zero temperature frontier between 
ordered and disordered phases, (A/J) c , is attracted to 
the percolation fixed point of the site diluted Ising model. 
Contrarily to results from a cluster variational analysis 
p|, the value of (A/J) c increases as p increases. Such 
a behavior is also predicted by standard mean-field ap- 
proximations. 
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